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ABSTRACT Coronaviruses contain a small envelope membrane protein with cation-selective ion channel activity mediated by its 
transmembrane domain (ETM). In a computational study, we proposed that ion channel activity can be explained by either of two 
similar ETM homopentameric transmembrane a-helical bundles, related by a ~50° rotation of the helices. Later, we tested this 
prediction, using site-specific infrared dichroism of a lysine-flanked isotopically labeled ETM peptide from the virus responsible for 
the severe acute respiratory syndrome, SARS, reconstituted in lipid bilayers. However, the data were consistent with the presence 
of a kink at the center of the ETM a-helix, and it did not fit completely either computational model. Herein, we have used native ETM, 
without flanking lysines, and show that the helix orientation is now consistent with one of the predicted models. ETM only produced 
one oligomeric form, pentamers, in the lipid-mimic detergent dodecylphosphocholine and in perfluorooctanoic acid. We thus report 
the correct backbone model for the pentameric a-helical bundle of ETM. The disruptive effects caused by terminal lysines probably 
highlight the conformational flexibility required during ion channel function. 
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Coronaviruses are enveloped viruses responsible for a 
variety of acute and chronic diseases (1). Their members have 
an envelope (E) protein, embedded in the membrane, which 
is critical for virion morphogenesis (2,3). These E proteins 
have also shown cation-selective ion channel activity, not 
only in the virus causative of severe acute respiratory syn- 
drome (SARS), but also in other members of the coronavirus 
genus (4,5). Ion channel activity is inhibited by the drug 
hexamethylene amiloride (HMA) in murine hepatitis virus 
and in human coronavirus 229E, and this inhibition was 
correlated with decreased viral replication in these viruses. 
Also, HMA did not have antiviral activity against a recom- 
binant murine hepatitis virus where E protein had been 
deleted (5), pointing to coronavirus E proteins as a pharma- 
cological target. 

In severe acute respiratory syndrome (SARS), the E 
protein has 76 amino acids, and contains one transmembrane 
a-helical domain (ETM), which is responsible for the 
observed ion channel activity (4,6) via the likely formation 
of a transmembrane a-helical bundle. In a recent computa- 
tional work (7), we suggested possible models for ETM 
oligomerization by exploring exhaustively the conforma- 
tional space of several ETM oligomers (8). This search was 
run in parallel with homologous sequences of ETM present 
in other coronaviruses, on the assumption that they share a 
common backbone structure (9). Conservative mutations that 
appear during evolution act as a filter, destabilizing nonna- 
tive models. At the same time, these mutations leave 
unaffected those models representing the native backbone 
structure. Under these evolutionary constraints, only two 
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models compatible with ion channel activity were found, 
both pentameric, referred to as models A and B (7). These 
two models had the same handedness and similar helix tilt, 
and were related by a ~50° rotation of their a-helices (7) 
(see Table 1, Models A and B). In 14 out of 15 ETM homol- 
ogous sequences, lowest-energy models clustered around 
model A (11), which was therefore likely to represent the 
native backbone structure of the ETM ion channel. 

To test our computational prediction, we performed 
recently a site-specific infrared dichroism study (12) of 
ETM in the SARS E sequence, reconstituted in dimyristoyl- 
phosphatidylcholine lipid bilayers. ETM was isotopically- 
labeled with '°C='*o (10), at nine native residues, except at 
F-23, where we used Ala because of the lack of labeled Phe. 
In addition, to facilitate purification (13), ETM was flanked 
by two lysines at each N- and C termini (K,-ETM-K;), as it 
has been shown for various a-TM oligomers that even four 
flanking lysines at each C- and N-terminus do not affect 
oligomerization during SDS-PAGE (14,15). Although the 
rotational orientation (12), w, of the centrally located labels 
was consistent with one of the models (B), labels located at 
the N- and C-ends of the TM were consistent with model A 
(10), suggesting a mild twist, or kink, in a region of the helix 
near residue F-23 (see Table 1, column K3-ETM-K, (A-23)). 

To test the possible effect of the mutation F-23A on ETM 
structure, we have studied herein the same peptide, K2-ETM- 
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K,, but using the native F-23 residue isotopically labeled 
with ?C=!8O, We note that the use of label F-23, instead of 
A-23, should not only affect the w value at F-23, but also the 
calculated value for the preceding and following residue. 
Indeed, the rotational orientation, w, is calculated by 
combining two consecutive residues (12), which are as- 
sumed to be separated by Aw = 100° for canonical a-helices. 
The infrared spectra obtained for K,-ETM-K, were similar to 
those obtained previously (10), and are not shown. 

Thus, from the measured dichroic ratios of the K,-ETM- 
K, sample, Rheiix and Rie (Shown in the Supplementary 
Material, Table 1, Data S1), we recalculated the orientational 
data for the labeled residues 22 and 24 (10). Table 1 (column 
K,-ETM-K, (F-23)) shows that the recalculated w values are 
still incompatible with either model A or B. 

Next, to test the effect of terminal lysine residues on ETM 
structure, we used the SARS ETM peptides shown in Fig. 1, 
where terminal lysines have been substituted to native ETM 
residues. These peptides were labeled independently at five 
native consecutive residues: L-21, A-22, F-23, V-24, and V-25. 

The data were analyzed as above (dichroisms are shown in 
the Supplementary Material, Table 1, Data S1), and the cal- 
culated w values are indicated in the last column of Table 
1 (ETM). These values are clearly consistent with model A; 
therefore, we conclude that the presence of terminal lysines in 
K,-ETM-K;j introduces a perturbation in the ETM structure, 
at least from residues 19 to 24. Consistent with this, solution 
NMR results of ETM in dodecylphosphocholine (DPC) 
micelles also indicate a kink in the central part of the a-helix 
in the peptide K-ETM-K3, but not in ETM (Fig. S1, Data S1). 

To confirm that ETM forms pentamers, we performed a 
sedimentation equilibrium study in DPC micelles. DPC is a 
mimic for dimyristoylphosphatidylcholine, as it has the same 
polar headgroups, and has been routinely used to determine 
oligomeric size of TM a-helices, e.g., (16). The sedimen- 


TABLE 1 

Ko-ETM-K> K>-ETM-K> 
Residue Model A Model B (A-23) (F-23) ETM 
V-17 31 80 AT = 5 
L-18 57 8 15: =: 8 
L-19 164 110 105 + 3 
F-20 —78 —132 
L-21 6 —30 —-44+4 4+4 
A-22 106 60 50 +7 220} 113 +8 
F-23 128 183 170 +7 -76+7 -109+7 
V-24 32 83 70+ 4 2444 -6+5 
V-25 80 22 94+ 5 
F-26 —173 130 
L-27 —64 —113 —48 + 4 
L-28 22 —15 52 


Rotational orientation values, w, for the predicted pentameric models A and 
B (7), experimental values from K2-ETM-K, using label A-23 (10), from 
K,-ETM-K, using native label F23 (this work), and native ETM without 
flanking lysines (this work). The helix tilt for the last two columns was 20— 
25° and is not shown. 
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ETGTLIVNSVLLELAFVVELLVILALLTALR 
ETGTLIVNSVLLFLAFVVFLLVTLAILTALR 


ETM ETGTLIVNSVLLFLAFVVFLLVTLATLTALR 
ETGTLIVNSVLLFLAFVVFLLVTILAILTALR 
ETGTLIVNSVLLFLAFVVFLLVTLATLTALR 

FIGURE 1 Synthetic peptide sequences of SARS ETM, with 


isotopically (‘SC='0) labeled residues indicated in bold, used 
for site-specific infrared dichroism. N- and C-terminal lysines 
used in a previous work (10) have been changed to the residues 
underlined, which correspond to native SARS ETM residues. 


tation data could be best fitted to a monomer-pentamer self- 
association model, with an association constant of 17.1 (Fig. 
$2, Data S1), although a higher order aggregate (n > 10) was 
also present (~ 10% of the material) at high concentrations of 
the peptide (not shown). 

We also performed electrophoresis of ETM in the mild 
detergent perfluorooctanoic acid (PFO). Although less 
suitable than zwitterionic DPC, this detergent has been 
successfully used to determine the oligomeric sizes of a-helical 
transmembrane peptides, e.g., in the tetrameric phospholem- 
man (17), or the pentameric small hydrophobic protein in 
the respiratory syncytial virus (18), which in SDS produced 
either monomers or multiple nonspecific oligomers, re- 
spectively. Fig. 2 shows that the only oligomer of ETM in 
PFO is pentameric (Jane 3). By contrast, in K2-ETM-K> 
(lane 2), the pentamer was not stable, and only produced 
dimers and trimers. Results in SDS were identical as those 
described previously (6,7,10) and are shown in Fig. S4 
(Data S1). 

These results suggest that the terminal lysines in K,-ETM- 
K, not only induce a kink in the ETM a-helix, but also may 
weaken interhelical interactions. These two effects, however, 
do not inactivate the ETM ion channel present in lipid 
bilayers, because K2-ETM-K, shows sodium ((6) and Fig. 
$5, Data S1) conductance activity, which was blocked by 
amantadine (6). ETM, however, displayed a different 
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FIGURE 2 Electrophoresis of ETM in PFO. Lane 1, molecular 
weight markers; lane 2, K2-ETM-K2; and lane 3, ETM without 
lysines(see Fig. 1). Arrows indicate monomers, dimers, and 
trimers (/ane 2), and monomers and pentamers (/ane 3). 
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conductance behavior relative to K,-ETM-K>, with an open 
probability that was proportional to the potential (Fig. S5, 
panel d, Data S1). By contrast, in K7-ETM-Kg, this depen- 
dency was lost. In addition, amantadine was less effective as 
an inhibitor for ETM than for K,-ETM-K; (Fig. S6, Data S1). 

Overall, our results show that flanking lysines induce a 
kink in ETM and weaken interhelical interactions. The 
precise origin of these effects is unknown, but it is likely 
related to the interaction of lysine side chains with the 
phosphocholine interface headgroups. These effects may 
unveil the presence of inherent conformational plasticity in 
ETM, possibly required during ion channel function, as 
suggested for other ion channels (19,20). Indeed, the central 
region of ETM, i.e., two turns from residue 19 to 26, contains 
three phenylalanines (bulky) and two valines (6-branched), 
which are known to destabilize a-helices (21). 

Our previous results (6) indicate that amantadine has more 
affinity for this lysine-induced conformational state than for 
the model represented by the native ETM, model A (10). 
This could explain the observed lack of protection of 
amantadine against SARS-induced cytopathic effects (22). A 
more detailed NMR study of the structure of this channel, 
complexed with the inhibitory drug HMA, is under way. 


SUPPLEMENTARY MATERIAL 


To view all of the supplemental files associated with this 
article, visit www.biophysj.org. 
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